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We present an updated measurement of the unitarity triangle angle 03 using a Dalitz plot analysis 
of the Kg-K + n~ decay of the neutral D meson produced in B ± —* D^K^ decays. The method 
exploits the interference between D° and D° to extract the angle 03, strong phase <5 and the ratio r of 
suppressed and allowed amplitudes. We apply this method to a 605 fb _1 data sample collected by the 
Belle experiment. The analysis uses two modes: B + — > DK + , and B + — ► D* K + with D* — ► Dty°, 
as well as the corresponding charge-conjugate modes. From a combined maximum likelihood fit to 
the two modes, we obtain 03 = 76° ^^(stat) ± 4°(syst) ± 9° (model). The statistical significance 
of CP violation (03 ^ 0) in our measurement is (1 — 5.5 x 10~ 4 ), or 3.5 standard deviations. These 
results are preliminary. 

PACS numbers: 12.15.Hh, 13.25.Hw, 14.40.Nd 



I. INTRODUCTION 

Determinations of the Cabibbo-Kobayashi-Maskawa 
(CKM) [1] matrix elements provide important checks 
on the consistency of the standard model and ways to 
search for new physics. The possibility of observing di- 
rect CP violation in B — > DK decays was first discussed 
by I. Bigi, A. Carter and A. Sanda [2]. Since then, var- 



ious methods using CP violation in B — ► DK decays 
have been proposed [3-6] to measure the unitarity trian- 
gle angle 03. Three body final states such as Kgir + ir~ 
[7, 8] have been suggested as promising modes for the 
extraction of 03 . In the Wolfenstein parameterization of 
the CKM matrix elements, the weak parts of the am- 
plitudes that contribute to the decay B + — * DK + are 
given by V* b V us - AX 3 (for the D°K+ final state) and 



3 



V* b V ca ~ AA 3 (p + irf) (for D°K+). The two amplitudes 
interfere as the Z?° and D° mesons decay into the same 
final state KgTr + Tr~. Assuming no CP asymmetry in 
neutral D decays, the amplitude of the neutral D decay 
from B ± — > DK^ as a function of Dalitz plot variables 
m 2 , = m 2 i and m 2 = m 2 a _ is 

M± = /(ml.mD + re^+^/Cm^mi), (1) 

where f(m 2 + , m 2 _) is the amplitude of the D° — > A!j7r + 7r~ 
decay, r is the ratio of the magnitudes of the two inter- 
fering amplitudes, and 5 is the strong phase difference 
between them. The D° -> K^tt+tt- decay amplitude / 
can be determined from a large sample of flavor-tagged 
D° — > A"g7r + 7r~ decays produced in continuum e + e~ an- 
nihilation. Once / is known, a simultaneous fit of B + 
and B~ data allows the contributions of r, 03 and (5 to be 
separated. The method has a two-fold ambiguity: (03 , <5) 
and (03 + 180°, 5 + 180°) solutions cannot be separated. 
We always choose the solution with < 03 < 180°. Ref- 
erences [7] and [9] give a more detailed description of the 
technique. 

The method described above can be applied to other 
modes as well as B + — > DK + decay and its charge- 
conjugate mode (charge conjugate states are implied 
throughout the paper). Excited states of neutral D 
and K mesons can also be used, although the values of 
8 and r can differ for these decays. Both BaBar and 
Belle collaborations have successfully applied this tech- 
nique to B ± — » £)(*)ii'(*) ± modes with D° decaying to 
KgTT + TT~ [9-f3]. In addition, the BaBar collaboration 
reported a measurement of 03 using the B^ — > DK ± 
mode with D° decaying to the 7r°7r + 7r _ final state [14]. 

In the current paper, we report a measurement of 03 
using the modes B+ -> DK+ and B+ -» D*K+ with 
D° — > KgTr + ir~ , based on a 605 fb _1 data sample col- 
lected by the Belle detector at the KEKB asymmetric 
e + e~ factory. These results are preliminary. 

II. EVENT SELECTION 

We use a data sample of 657 x 10 6 BB pairs, collected 
by the Belle detector. The decay chains B + — * DK + 
and B+ -* D*K+ with D* -> Dtt are selected for the 
analysis. The neutral D meson is reconstructed in the 
KgTr + ir~ final state in all cases. We also select decays of 
D*~ — > D°7r~ produced via the e + e~ —> cc continuum 
process as a high-statistics sample to determine the D° — ► 
Kgir + ir~ decay amplitude. 

The Belle detector is described in detail elsewhere 
[15, 16]. It is a large-solid-angle magnetic spectrometer 
consisting of a silicon vertex detector (SVD), a 50-layer 
central drift chamber (CDC) for charged particle tracking 
and specific ionization measurement (dE/dx), an array 
of aerogel threshold Chercnkov counters (ACC), time- 
of- flight scintillation counters (TOF), and an array of 
CsI(Tl) crystals for electromagnetic calorimctry (ECL) 



located inside a superconducting solenoid coil that pro- 
vides a 1.5 T magnetic field. An iron flux return located 
outside the coil is instrumented to detect Kl mesons and 
identify muons (KLM). 

Charged tracks are required to satisfy criteria based 
on the quality of the track fit and the distance from the 
interaction point. We require each track to have a trans- 
verse momentum greater than 100 MeV/c. Separation 
of kaons and pions is accomplished by combining the re- 
sponses of the ACC and the TOF with the dE/dx mea- 
surement from the CDC. Photon candidates are required 
to have ECL energy greater than 30 MeV. Neutral pion 
candidates arc formed from pairs of photons with invari- 
ant masses in the range 120 to 150 MeV/c 2 . Neutral 
kaons are reconstructed from pairs of oppositely charged 
tracks with an invariant mass M nv within 7 MeV/c 2 of 
the nominal Kg mass and the reconstructed vertex dis- 
tance from the interaction point greater than 1 mm. 

To determine the D° -> K s ti + h^ decay amplitude we 
use D*^ mesons produced via the e + e~ — > cc continuum 
process. The flavor of the neutral D meson is tagged by 
the charge of the slow pion (which we denote as n s ) in 
the decay D*~ — > D°ir~ . The slow pion track is required 
to originate from the D° decay vertex to improve the 
momentum and angular resolution of the n s . To select 
neutral D candidates we require the invariant mass of 
the A^tt+tt- system to be within 11 MeV/c 2 of the D° 
mass. To select events originating from a D*~ decay we 
impose a requirement on the difference AM of the in- 
variant masses of the D*~ and the neutral D candidates: 
144.9 MeV/c 2 < AM < 145.9 MeV/c 2 . The suppres- 
sion of the combinatorial background from BB events is 
achieved by requiring the D*~ momentum in the center- 
of-mass (CM) frame to be greater than 2.7 GeV/c. The 
number of events in the signal region is 290.9 x 10 3 ; the 
background fraction is 1.0%. 

The selection of B candidates is based on the CM 
energy difference AE = ^ Ei — -Ebcam and the beam- 
constrained B meson mass Mb c = v^beam ~ (Eft) 2 ! 
where £?beam is the CM beam energy, and Ei and pi are 
the CM energies and momenta of the B candidate decay 
products. We also impose a requirement on the invariant 
mass of the neutral D candidate: \M K o^+^- — M D a\ < 
11 MeV/c 2 . 

To suppress background from e + e~ — > qq (q = 
u, d, s, c) continuum events, we calculate two variables 
which characterize the event shape. One is the cosine of 
the thrust angle cos#thr, where 6* t hr is the angle between 
the thrust axis of the B candidate daughters and that of 
the rest of the event. The other is a Fisher discriminant 
T composed of 11 parameters [17]: the production angle 
of the B candidate, the angle of the B thrust axis rela- 
tive to the beam axis and nine parameters representing 
the momentum flow in the event relative to the B thrust 
axis in the CM frame. At the first stage of the analy- 
sis, when the (Mbc AE) distribution is fitted in order to 
obtain the fractions of the background components, we 
require | cos(9 t hr| < 0.8 and T > —0.7. In the Dalitz plot 
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FIG. 1: AE and M bc distributions for the B + -> DA"+ (top) 
and B + — * D* K + (bottom) event samples. Points with er- 
ror bars are the data, and the histogram is the result of a 
MC simulation according to the fit result. The AE (Mbc) 
distributions are shown here with a signal-region selection of 
M bc > 5.27 GeV/c 2 (\AE\ < 30 MeV) applied; this fit is 
performed on the full region. 



fit, we do not reject events based on these variables (as 
in the previous analysis [9]), but rather use them in the 
likelihood function to better separate signal and back- 
ground events. This leads to a 7-8% improvement in the 
expected statistical error. 

The AE and M hc distributions for B + -> DK + and 
B + — » D*K + candidates are shown in Fig. 1. For the se- 
lected events a two-dimensional unbinncd maximum like- 
lihood fit in the variables Mb c and AE is performed, with 
the fractions of continuum, BB and i? ± — > D^*V ± back- 
grounds as free parameters, and their distributions fixed 
from generic MC simulation. The resulting signal and 
background density functions are used in the Dalitz plot 
fit to obtain the event-by-event signal to background ra- 
tio. The number of events in the signal box (Mt, c > 5.27 
GeV/c 2 , \AE\ < 30 MeV, |cos0 t hr| < 0.8, T > -0.7) is 
756. The (Mbc AiS) fit yields a continuum background 
fraction of (17.9 ± 0.7)%, BB background fraction of 
(7.3 ± 0.5)%, and a B ± — > Dir^ background fraction 
of (4.3 ± 0.3)% in the signal box. 

To select B + — > D*K + events, in addition to the re- 
quirements described above, we require that the mass 
difference AM of neutral D* and D candidates satis- 
fies 140 MeV/c 2 < AM < 144 MeV/c 2 . The number of 
events in the signal box is 149. The continuum back- 
ground fraction is (5.7 ±0.7)%, the BB background frac- 
tion is (7.6 ± 1.9)%, and B ± — > L»*7r ± background frac- 
tion is (7.0 ± 1.3)%. 



FIG. 2: Dalitz distributions of D° J\& n + ir decays from 
selected B ± -> DA"* (top) and B ± -> D*K ± (bottom) can- 
didates, shown separately for B~ (left) and B + (right) tags. 



The Dalitz distributions of D° — > if!j7r + 7r _ decay in 
the signal box of B ± — > Dif ± and B ± — > D*K ± pro- 
cesses are shown in Fig. 2. 



III. DETERMINATION OF THE D° -> K^tt+tt- 
DECAY AMPLITUDE 

As in our previous analysis [9], the D° — > Kgir + ir~ 
decay amplitude is represented using the isobar model. 
The list of resonances is also the same, the only dif- 
ference being the free parameters (mass and width) of 
the if* (892) * and p(770) states. A modified amplitude, 
where the scalar irir component is described using the 
K-matrix approach [18], is used in the estimation of the 
systematic error. 

The amplitude / for the D° — > KgTT + TT~ decay is de- 
scribed by a coherent sum of N two-body decay ampli- 
tudes and one non-resonant decay amplitude, 

jv 

/ (m\ ,m 2 _) = Y, a 3 Aj (m\ , m 2 _ ) + a NR e« NR , (2) 
j'=i 

where Aj{m\,rn_) is the matrix clement, aj and £j 
are the amplitude and phase of the matrix element, 
respectively, of the j-th resonance, and onr and £nr 
are the amplitude and phase of the non-resonant com- 
ponent. The description of the matrix elements fol- 
lows Rcf. [19]. We use a set of 18 two-body am- 
plitudes. These include five Cabibbo-allowed am- 
plitudes: if*(892)+7r-, if*(1410)+7r-, if*(1430)+7r-, 
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jq(1430)+7r- and ^*(1680) + tt-; their doubly Cabibbo- 
suppressed partners; and eight amplitudes with Kg and 
a 7T7T resonance: K° s p, K^lv, A^/ o (980), ifg/ 2 (1270), 
ifg/o(1370), A^p(1450), ifgo-i and K° s a 2 . 

We use an unbinned maximum likelihood technique to 
fit the Dalitz plot distribution to the model described by 
Eq. 2 with efficiency variation, background contributions 
and finite momentum resolution taken into account. The 
free parameters of the minimization are the amplitudes 
dj and phases £j of the resonances, the amplitude a^R 
and phase £nr of the non-resonant component and the 
masses and widths of the <j\ and a 2 scalars. We also 
allow the masses and widths of the K* (892) + and p(770) 
states to float. 

The procedures for determining the background den- 
sity, the efficiency, and the resolution are the same as in 
the previous analyses [9, 11]. The background density for 
D — > Kgir + ir~ events is extracted from AAI sidebands. 
The shape of the efficiency over the Dalitz plot, as well 
as the invariant mass resolution, is extracted from the 
signal Monte-Carlo (MC) simulation. 

The fit results are given in Table I. The parameters 
obtained for the o x resonance (M ai = 522 ± 6 MeV/c 2 , 
r CTl = 453 ± 10 MeV/c 2 ) are similar to those observed by 
other experiments [20, 21]. The second scalar term (T2 is 
introduced to account for structure observed at to 2 ,,. ~ 
1.1 GeV 2 /c 4 : the fit finds a small but significant contribu- 
tion with M a2 = 1033 ±7 MeV/c 2 , T a2 = 88 ±7 MeV/c 2 . 
Allowing the parameters of the dominant if *(892) + and 
p(770) resonances to float results in a significant im- 
provement in the fit quality. We obtain M(K* (892)) = 
893.7 ± 0.1 MeV/c 2 , T(K*(892)) = 48.4 ± 0.2 MeV/c 2 , 
M(p) = 771.7±0.7MeV/c 2 , T(p) = 136.0 ± 1.3 MeV/c 2 . 

The x 2 test finds a value of x 2 / n df = 2.35 for 1065 
degrees of freedom (ndf), which is large. We find that 
the main features of the Dalitz plot are well-reproduced, 
with some significant but numerically small discrepancies 
at peaks and dips of the distribution. In our final results 
we include a conservative contribution to the systematic 
error due to uncertainties in the D° decay model. 

IV. DALITZ PLOT ANALYSIS OF B+ -> D^ ) K + 
DECAYS 

As in our previous analysis [9] and in analyses car- 
ried out by the BaBar collaboration [12, 13], we fit the 
Dalitz distributions of the B + and B~ samples sepa- 
rately, using Cartesian parameters x± = r± cos(±03 + S) 
and y± = r± sin(±03 + 5), where the indices "+" and 
"— " correspond to B + and B~ decays, respectively. In 
this approach the amplitude ratios (r+ and r_) are not 
constrained to be equal for the B + and B~ samples. Con- 
fidence intervals in r, 4>s and 5 are then obtained from 
the (x±,y±) using a frequentist technique. The advan- 
tage of this approach is low bias and simple distributions 
of the fitted parameters, at the price of fitting in a space 
with higher dimensionality (2;+, y+, X-, yJ) than that of 



TABLE I: Fit results for D° -> iTg7r + 7r~ decay. Errors are 
statistical only. 



Intermediate state 


Amplitude 


Phase (°) 


Ks(T\ 


1.56 ±0.06 


214 ±3 


Ksp° 


1.0 (fixed) 


(fixed) 


K s oj 


0.0343 ± 0.0008 


112.0 ± 1.3 


K s /o(980) 


0.385 ± 0.006 


207.3 ± 2.3 


Kscr 2 


0.20 ± 0.02 


212 ± 12 


A 5 / 2 (1270) 


1.44 ±0.04 


342.9 ± 1.7 


if s /o(1370) 


1.56 ±0.12 


110 ±4 


K s p°(1450) 


0.49 ± 0.08 


64 ±11 


A'*(892)+7T- 


1.638 ±0.010 


133.2 ±0.4 


A*(892)"7r+ 


0.149 ± 0.004 


325.4 ± 1.3 


A*(1410) + tt- 


0.65 ± 0.05 


120 ±4 


A*(1410)~7T+ 


0.42 ± 0.04 


253 ±5 


Ao(1430)+7r- 


2.21 ±0.04 


358.9 ±1.1 


^(1430)-tt + 


0.36 ± 0.03 


87 ±4 


A 2 *(1430)+7r- 


0.89 ± 0.03 


314.8 ± 1.1 


A 2 *(1430)~tt + 


0.23 ± 0.02 


275 ±6 


A*(1680)+7r- 


0.88 ± 0.27 


82 ± 17 


A*(1680)-7T+ 


2.1 ±0.2 


130 ±6 


non-resonant 


2.7 ±0.3 


160 ±5 



the physical parameters (r,(j>3,6); see Section V. 

The fit to a single Dalitz distribution with free param- 
eters x and y is performed by the unbinned maximum 
likelihood technique, using variables m 2 + , to 2 , Mb c , AE, 
cos 6> t hr and T\ only the first four variables were used in 
the previous analysis [9]. We subdivide the background 
distribution into four components: e + e~ — > qq (where 
q = u,d,s), charm, BB (except for B ± — ► D^ix^) and 
£? ± — ► D^'if^- background. The distributions of each of 
these components are assumed to be factorized as prod- 
ucts of a Dalitz plot distribution (m 2 ^, to 2 .), and distribu- 
tions in (M bc ,AE), and (cos t hr > F) ■ The shapes of these 
distributions are extracted from MC simulation. The six- 
dimensional PDF used for the fit is thus expressed as 

p = Y J P l (m 2 +, m 2 _ )pt (M hc , AE) Pi (cos 6> thr , F), (3) 

i 

where the index i runs over all background contributions 
and signal. Possible deviations from the factorization 
assumption and disagreements between MC and experi- 
mental background densities are treated in the systematic 
error. The efficiency variation as a function of the Dalitz 
plot variables is obtained from signal MC simulation and 
is taken into account in the likelihood function. 

To test the consistency of the fit, the same procedure as 
used for B + — » D^K + signal was applied to the B + — > 
£>(*'7r + control samples. The results are consistent with 
the expected value r ~ 0.01 for the amplitude ratio. 
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TABLE II: Results of the signal fits in parameters (x, y). The 
first error is statistical, the second is experimental systematic 
error. Model uncertainty is not included. 



Parameter 


B+ -> DK+ 


B+ -> D*K+ 


X- 


+0.105 ±0.047 ±0.011 


±0.024 ±0.140 ±0.018 


y- 


±0.177 ±0.060 ±0.018 


-0.243 ±0.137 ±0.022 


x+ 


-0.107 ±0.043 ±0.011 


±0.133 ±0.083 ±0.018 


y+ 


-0.067 ±0.059 ±0.018 


±0.130 ±0.120 ±0.022 




-0.2 0.2 -0.2 0.2 



FIG. 3: Results of signal fits with free parameters x = rcosO 
and y = rsin(9 for B+ -> DK+ (left), and B+ -> D*K+ 
(right) samples, separately for B~ and B + data. Contours 
indicate one and two standard deviation regions obtained in 
the maximum likelihood fit. 




50 100150 200 250 300 350 50 100150 200 250 300 350 

o 3 (degrees) o 3 (degrees) 



FIG. 4: Projections of confidence regions for the B + — * DK + 
(top) and B + — * D*K + (bottom) mode onto the (r,fa) and 
(fa, 5) planes. Contours indicate projections of one, two and 
three standard deviation regions. 



VI. ESTIMATION OF SYSTEMATIC ERROR 



The results of the separate B + and B~ data fits are 
shown in Fig. 3. The values of the fit parameters x± and 
y± are listed in Table II. 



V. EVALUATION OF THE STATISTICAL 
ERRORS 

We use a frequentist technique to evaluate the statis- 
tical significance of the measurements. The procedure 
is identical to that in our previous analysis [9], This 
method requires knowledge of the probability density 
function (PDF) of the reconstructed parameters x and 
y as a function of the true parameters x and y. To ob- 
tain this PDF, we employ a "toy" MC technique that 
uses a simplified MC simulation of the experiment which 
incorporates the maximum likelihood fit with the same 
efficiencies, resolution and backgrounds as used in the fit 
to the experimental data. 

Figure 4 shows the projections of the three-dimensional 
confidence regions onto the (r, fa) and (fa, 6) planes for 
B ± -> DK± and B ± -> D*K ± modes. We show the 
20%, 74% and 97% confidence level regions, which cor- 
respond to one, two, and three standard deviations for a 
three-dimensional Gaussian distribution. The values of 
the parameters r, fa and (5 obtained for B — » DK^ 
and B ± — > D*K ± modes separately are presented in Ta- 
ble III. 



Experimental systematic errors come from the uncer- 
tainty in the knowledge of the distributions used in the 
fit (i.e. Dalitz plot distributions of the background com- 
ponents, and the (M bc ,AE) and (cos# thr , T) distribu- 
tions of the backgrounds and signal), fractions of dif- 
ferent background components, and the distribution of 
the efficiency across the Dalitz plot phase space. Un- 
certainties due to background shapes are estimated by 
using alternative distributions in the fit (extracted from 
experimental data where possible). Uncertainties due to 
the background fractions are obtained by varying each 
fraction within its error. The procedure for estimating 
the uncertainty due to the detection efficiency is modi- 
fied compared to the previous analysis [9] : we use an al- 
ternative efficiency shape obtained by toy MC from the 
parameterized track finding efficiency (obtained from ex- 
perimental data) as a function of transverse momentum 
and azimuthal angle 6. 

Compared to our previous analysis [9], an additional 
source of systematic error exists due to the use of cos # t hr 
and J- variables in the fit. However, the use of these 
variables increases the effective signal-to-background ra- 
tio, so the total systematic error is comparable. 

Systematic errors in the physical parameters r, fa and 
8 are calculated from the systematic errors on the fitted 
parameters (x, y). Values {x,y) are generated accord- 
ing to Gaussian distributions with standard deviations 
equal to the corresponding total systematic errors, then 
parameters r, fa and 5 are obtained for each (x, y) set, 
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TABLE III: CP fit results. The first error is statistical, the second is experimental systematic, and the third is the model 
uncertainty. 



Parameter 


B + -> DK+ mode 


B + -> D*K + mode 


03 

r 
8 


80.8° tllil ± 5.0° ± 8.7° 
0.16lio.o3g ±0.011 ±0.049 
137.4° t^-o» ±4 .o° ±2 2.9° 


63.8° j° ± 4.7° ± 8.7° 

0.208to:°83 ± °- 015 ± °- 049 
342.0° t^jo ± 3.7° ± 22.9° 



and the root-mean-square deviations (RMS) of the re- 
sulting values are calculated. We perform this procedure 
in two ways: without correlation of (x,y) biases for B + 
and B~ , and with 100% correlation between them. The 
largest RMS of the two options serves as the systematic 
error. The systematic errors in x, y variables are shown 
in Table II. 

The model used for the D° — > KgTr + ir~ decay ampli- 
tude is one of the main sources of error for our analy- 
sis: we list this contribution separately. For the current 
measurement we use the same estimates of the model un- 
certainty as in our previous analysis [9]. In addition, we 
perform a fit of the D° — > Kg ir + ir~ amplitude using the 
K-matrix formalism [18] to describe the irir S-w&ve con- 
tribution. The maximum difference between the baseline 
quasi two-body amplitude and the K-matrix amplitude 
in the 03 fit is ~ 2°. However, since the K-matrix de- 
scribes only part of the amplitude, we still use A03 = 9° 
as the estimate of the model uncertainty. 

VII. COMBINED 3 MEASUREMENT 

The two event samples, B + — > DK + and B + — » 
D*K + are combined in order to improve the sensitiv- 
ity to 03. The confidence levels for the combination of 
two modes are obtained using the frequentist technique 
as for the single mode, with the PDF of the two measure- 
ments being the product of the probability densities for 
the individual modes. Confidence intervals for the com- 
bined measurement together with systematic and model 
errors are shown in Table IV. The statistical confidence 
level of CP violation is (1 — 5.5 x 10~ 4 ), or 3.5 standard 
deviations. 



VIII. CONCLUSION 

We report the results of a measurement of the unitar- 
ity triangle angle 03, using a method based on Dalitz 



plot analysis of D° — > Alj7r + 7r _ decay in the process 
B + — ► D^*'K + . An updated measurement of 03 us- 
ing this technique was performed based on 605 fb _1 of 
data collected by the Belle detector: 70% larger than the 
previous sample. The statistical sensitivity of the mea- 
surement has also been improved by modifications to the 
event selection and fit procedure. 

From the combination of B + — ► DK + and B + — » 
D*K + modes, we obtain the value 03 = 76° + 13 o (stat) ± 
4°(syst) ± 9° (model); of the two possible solutions we 
choose the one with < 03 < 180°. We also obtain val- 
ues of the amplitude ratios tdk = 0-16 ± 0.04(stat) ± 
O.Ol(syst) ± 0.05(model) and r D , K = 0.21 ± O.OS(stat) ± 
0.02(syst) ± 0.05(model). The statistical significance of 
CP violation for the combined measurement is (1 — 5.5 x 
10~ 4 ), or 3.5 standard deviations. These results are pre- 
liminary. 

The statistical precision of the 03 measurement is al- 
ready comparable to the estimated model uncertainty. 
However, it is possible to eliminate this uncertainty us- 
ing decays of ^(3770) -> D°D° [22, 23]. 



Acknowledgments 

We thank the KEKB group for excellent operation of 
the accelerator, the KEK cryogenics group for efficient 
solenoid operations, and the KEK computer group and 
the Nil for valuable computing and SINET3 network sup- 
port. We acknowledge support from MEXT and JSPS 
(Japan); ARC and DEST (Australia); NSFC and KIP 
of CAS (China); DST (India); MOEHRD, KOSEF and 
KRF (Korea); KBN (Poland); MES and RFAAE (Rus- 
sia); ARRS (Slovenia); SNSF (Switzerland); NSC and 
MOE (Taiwan); and DOE (USA). 



[1] M. Kobayashi and T. Maskawa, Prog. Theor. Phys. 49, 
652 (1973); N. Cabibbo, Phys. Rev. Lett. 10, 531 (1963). 

[2] I. I. Bigi and A. I. Sanda, Phys. Lett. B211, 213 (1988); 
A. B. Carter and A. I. Sanda, Phys. Rev. Lett 45, 952 
(1980). 



[3] M. Gronau and D. London, Phys. Lett. B253, 483 
(1991); M. Gronau and D. Wyler, Phys. Lett. B265, 172 
(1991). 

[4] I. Dunietz, Phys. Lett. B270, 75 (1991). 

[5] D. Atwood, G. Eilam, M. Gronau and A. Soni, Phys. 



8 



TABLE IV: Results of the combination of B + -> DK + and B + -> ZT.ft~ + modes. 



Parameter 


la interval 


2a interval 


Systematic error 


Model uncertainty 




yeo +12° 
fU _ 13 o 


49° < 4> B < 99° 


4° 


9° 


fDJT 


0.16 ±0.04 


0.08 < r DK < 0.24 


0.01 


0.05 


VD*K 


0.21 ±0.08 


0.05 < r D -K < 0.39 


0.02 


0.05 


Sdk 


136° 


100° < 8 DK < 163° 


4° 


23° 


Sd*k 


343° ±™l 


293° < 8dk < 389° 


4° 


23° 



Lett. B341, 372 (1995). 
[6] D. Atwood, I. Dunietz and A. Soni, Phys. Rev. Lett. 78, 
3257 (1997); D. Atwood, I. Dunietz and A. Soni, Phys. 
Rev. D 63, 036005 (2001). 
[7] A. Giri, Yu. Grossman, A. Soffer, J. Zupan, Phys. Rev. 

D 68, 054018 (2003). 
[8] A. Bondar. Proceedings of BINP Special Analysis Meet- 
ing on Dalitz Analysis, 24-26 Sep. 2002, unpublished. 
[9] Belle Collaboration, A. Poluektov et al., Phys. Rev. D 
73, 112009 (2006). 
[10] Belle Collaboration, K. Abe et al., hep-ex/0308043. 
[11] Belle Collaboration, A. Poluektov et al, Phys. Rev. D 

70, 072003 (2004). 
[12] BaBar Collaboration, B. Aubert et al., Phys. Rev. Lett. 

95, 121802 (2005). 
[13] BaBar Collaboration, B. Aubert et al., hep-ex/0507101. 
[14] BaBar Collaboration, B. Aubert et al., Phys. Rev. Lett. 
99, 251801 (2007). 



[15] Belle Collaboration, A. Abashian et al., Nucl. Instr. and 

Meth. A 479, 117 (2002). 
[16] Y. Ushiroda (Belle SVD2 Group), Nucl. Instr. and Meth. 

A 511, 6 (2003). 
[17] CLEO Collaboration, D. M. Asner et al, Phys. Rev. D 

53, 1039 (1996). 
[18] V. V. Anisovich and A. V. Sarantsev, Eur. Phys. J. A 

16, 229 (2003). 

[19] CLEO Collaboration, S. Kopp et al., Phys. Rev. D 63, 
092001 (2001). 

[20] CLEO Collaboration, H. Muramatsu et al., Phys. Rev. 
Lett. 89, 251802 (2002), Erratum- ibid: 90, 059901 
(2003). 

[21] E791 Collaboration, E. M. Aitala et al, Phys. Rev. Lett. 
86, 770 (2001). 

[22] A. Bondar, A. Poluektov, Eur. Phys. J. C 47, 347 (2006). 
[23] A. Bondar, A. Poluektov, arXiv:0801.0840vl [hep-ex]. 



